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The results of a theoretical investigation on the stopping power of ions moving in a disordered 
, two-dimensional degenerate electron gas are presented. The stopping power for an ion is calculated 

employing linear response theory using the dielectric function approach. The disorder, which leads 
to a damping of plasmons and quasiparticles in the electron gas, is taken into account through a 
I *~ \ relaxation time approximation in the linear response function. The stopping power for an ion is 

calculated in both the low- and high-velocity limits. In order to highlight the effects of damping we 
present a comparison of our analytical and numerical results, in the case of point-like ions, obtained 
for a non-zero damping with those for a vanishing damping. It is shown that the equipartition 
sum rule first formulated by Lindhard and Winther for three-dimensional degenerate electron gas 
does not necessarily hold in two-dimensions. We have generalized this rule introducing an effective 
dielectric function. In addition some new results for two-dimensional interacting electron gas have 
I , been obtained. In this case the exchange-correlation interactions of electrons are considered via 

pH ■ local-field-corrected dielectric function. 

PACS numbers: 52.40.Mj, 52.25.Mq, 73.50.Mx, 52.27.Gr 

11: 

"55 ! I- INTRODUCTION 

, There is an ongoing interest in the theory of interaction of swift charged projectiles with condensed matter. Although 
most theoretical works have reported on the energy loss of ions in a target medium which is modelled as a three- 
dimensional (3D) electron gas, the two-dimensional (2D) case has not yet received as much attention as the 3D case. 
A 2D electron system is now experimentally realizable in a laboratory. In the last three decades or so many interesting 
and intriguing properties of a 2D electron gas have been explored. For a recent update on some of these developments 
we refer to Refs. [3, 0. A widely used 2D electron system is realized at the interface between GaAs and Gai-^Al^As, 
^» j and in the interface metal-oxide-semiconductors (MOS). The interaction of charged particles with an electron gas is 
^vq . an important probe of many-body interactions in the target electron medium. It is known that many-body properties 
of an electron gas vary in notable aspects with spatial dimensions. It is therefore of interest to make a detailed 
study of interaction of charged particles with a 2D electron gas. This theoretical study is also of relevance to device 
' applications e.g. in using ion implantation in devices which involve 2D electron systems. 

In general, interaction of charged projectiles with condensed matter can be studied by means of the stopping power 
(SP) of the target medium. The SP accounts for the energy loss by an external charged projectile as it passes through 
and interacts with matter. And the SP of a medium can be used to construct diagnostic tools for studying this kind 
of physical systems. There have been several theoretical approaches to the energy loss and SP for 3D systems, and 
some of these approaches have been applied also to a 2D electron gas. Amongprevious theoretical works on a 2D 
electrongas some are based on the linear response dielectric function method H, [f| 0, [f| and quantum scattering 
theory Q [Hi [HI [HJ . Further works have dealt with some nonlinear screening effects through a quadratic response 
approach within the random-phase approximation [l3j], the employment of density functional theory flij and in a 
method based on frequency moments of the energy loss function Q . 

In this paper we shall consider fast charged projectiles and hence a linear response theory to calculate energy loss 
is expected to be adequate. Previously, within this approach, Bret and Deutsch calculated the SP of an ion 0, Q 
and a dicluster [5j] in a 2D electron gas for any degeneracy. Their results show some interesting differences with 
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the corresponding results for a 3D case. Of special interest is their finding that the leading term of the asymptotic 
expansion of the SP in a high- velocity limit decreases as l/v, where v is the projectile velocity, which differs from the 
well-known form predicted by the Bethe-Bloch formula 0,13,0] m the 3D case. The calculations in Refs. 0,13,0] are 
based on the random-phase approximation (RPA) which works well if electron-electron interaction can be neglected. 
Now, in 2D systems, electron density can be varied. For moderate values of electron density e.g. in semiconductors 
electron-electron interaction may not be negligible and going beyond RPA is desirable. 

Our objective is to consider two physically motivated aspects of a 2D electron gas in the context of energy loss. 
For the first part of our study we consider a disordered electron gas which contains impurities. The effect of these 
impurities is included through a phenomenological relaxation time for electrons due to scattering by impurities. For 
this disordered system we use a linear response dielectric function in RPA and in a number-conserving relaxation time 
approximation (RTA), which was first considered by Mermin [l8[ and then by Das for a 3D electron gas. This RTA 
formulation has not yet been extended beyond RPA. The effect of disorder which leads to a damping of excitations 
enters the RPA dielectric function, for a given electron-impurity collision frequency, through £rpa(&, w + vy), where 
7 is used as a model parameter. For a degenerate electron gas (DEG) and for a given electron density, the damping 
parameter can be assumed to be a constant to a good approximation. The disorder-inclusive dielectric function, 
with the collision frequency as a free parameter, allows some physical insight and useful numerical estimates of 
the influence of disorder on energy loss in a DEG. In 3D the predicted effect is a shorter life time with a smaller 
propagation wavelength of plasmons resulting considerable modifications of the SP (see, e.g., Refs. [20l l2ll I22I [23l. l2"j] 
and references therein). For the stopping of a single ion, the broadening of the plasmon peak with increasing 7 shifts 
the threshold for energy loss by plasmon excitation towards lower projectile velocities. It now becomes possible for 
low-velocity projectile ions to excite plasmons (in addition to single-particle excitations ). T his increases the SP of 3D 
electron gas at low projectile velocities, compared to the disorder- free RPA result (22I l23l |24| . The situation with a 
2D electron gas will be discussed in detail in the following sections. 

The second objective of our study is to investigate the influence of exchange-correlation interaction (i.e. beyond 
RPA) in an electron gas on the SP. For a 3D electron system it has been shown [25|, [2(| that the SP in low and 
intermediate velocity regimes shows a definite increase due to this interaction. A similar result has been reported 
for a 2D system 0, 0] • However let us note that if an asymptotic expansion of the SP in a high- velocity regime is 
considered then it has been shown previously that the first term in this expansion is unaffected by electron-electron 
interaction. In this paper we calculate the next non-vanishing term of this asymptotic expansion and show that 
it behaves as B(r s )v~~ i ln[A(r s )v], where exchange-correlation interactions are involved in A(r a ) and B(r s ). These 
functions depend on the target density through Wigner-Seitz density parameter r s — {imQafy^ 1 / 2 , where n and ao 
are electron gas density and Bohr radius, respectively. The details are presented in Sec. IIII1 

The plan of the paper is as follows. In Sec. Ill Al we derive analytical expressions for the disorder-inclusive dielectric 
function (DF) for a 2D degenerate electron gas (DEG). We would like to mention that an alternative but equivalent 
derivation is presented in Appendix [A] The latter derivation contains certain attractive features. Through this 
alternative formulation we consider a small-fc,cj approximation for the DF, and this approximate result is used in 
Sec. Ill Al In a small-fc,c<j approximation the plasmon dispersion for a disordered two-dimensional DEG exhibits a 
constraint not present in 3D. This behavior has been previously discussed in the literature [27|, \2B, l29( . We revisit this 
approximation through our formulation in Appendix [Aj The exact plasmon dispersion relations for an interacting 
DEG (including exchange-correlation effects) are derived in Sec. Ill Bl bv employing local- field corrections to the RPA 
dielectric function. In Sec. IIIII we briefly outline the general linear response function formalism of the 2D stopping 
power of a point-like ion. After dealing with the excitation equipartition in Sec. IIII Al we develop, in Sees. IITTB1 and 
IIII CI some analytical techniques to calculate the SP of an ion in low- and high- velocity regimes. The two particular 
cases studied in these sections are (i) low-velocity limit of the SP for an ion moving in a disordered DEG, and (ii) 
high- velocity limit for a strongly interacting DEG. Sec. IIVI contains systematic numerical calculations for the SP. The 
results are summarized in Sec. |V] which also includes discussion and outlook. Appendix [A] to which we draw the 
reader's attention presents the above-mentioned alternative derivation of the DF for the disorder- inclusive case in 
RPA, which is also valid in the complex cj-plane. In Appendix [B] we provide some technical details for an evaluation 
of the asymptotic SP. 

II. DIELECTRIC FUNCTION AND DISPERSION RELATIONS FOR 2D ELECTRON GAS 

In the linear response theory, the stopping power (SP) of an external projectile moving in a medium is related to 
the dielectric function e(fc,w) of the medium. Both the single-particle and collective excitations (i.e. the plasmons) 
contribute to the SP and these contributions are contained in e{k,u)) (see, e.g., Eq. (f2"tJ|) below). In our study the 
two-dimensional (2D) target medium is assumed to be disordered due to impurities etc. We shall incorporate effects 
of disorder in e(k, to) in a somewhat phenomenological manner. This is to include disorder through a relaxation time 
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r such that the particle number is conserved. For a three-dimensional (3D) medium this was done first by Mermin 
fl8| and then by Das fli^| in the random phase approximation (RPA) and in relaxation time approximation (RTA). 
We refer the reader to [HI, [l^l for details of this formalism. For r — ► oo, this linear response function e(k,u>, 1/r) 
reduces to the Lindhard dielectric function [Tf|[l7|- The dielectric function e(k, u>) is understood to contain 7 (= 1/r) 
as a damping parameter due to disorder. The form of e{k,uj, 1/r) is to be specified shortly for a 2D electron gas. 

It is convenient to introduce the dimensionless Lindhard variables z — fc/2fcf, u — co/kvp, where vf and Uf 
(= (27T?io) 1 ^ 2 ) are, respectively, the Fermi velocity and wave number of the target electrons. Also we introduce the 
density parameters \ 2 = \ j a = l/kpaa = r s /V2- In our calculations x an d ct (or r s ) serve as a measure of electron 
density. (Note that the density parameter \ introduced above differs from usual definition by a factor 7r see, e.g., 
Refs. I, i, Si). 



A. Disordered electron gas: RPA 

Let us now specify the disorder-inclusive dielectric function for 2D zero-temperature (degenerate) electron gas 
(DEG). This has been done previously in Refs. [U [H, [29| employing small-fc,o; approximation. Here within RPA 
and RTA we derive the disorder-inclusive dielectric function (DF) without further approximations on the energy- 
momentum spectrum i.e. on u> and k. As pointed out in Ref. [28| the physical arguments for deriving number- 
conserving DF by Mermin [l8[ and Das [19( in 3D are independent of dimensionality. Therefore with the notations 
introduced in the preceding paragraph, the DF for 2D DEG reads 



s(z, u) = 1 + 



(zu + iT) [£rpa(z, u, T) - 1] 



iT [e RPA (z, u, T) - 1] / [e KPA (z, 0) - 1] ' 



(1) 



where T = hj/AEp, Ep being the Fermi energy = h 2 k F /2m with m as the effective mass. The quantity 7 (or T) is a 
measure of damping of excitations in the disordered electron gas. £rpa(z, u, T) = £RpA(fc, w + vy) is the longitudinal 
dielectric function of DEG in the RPA derived in 2D by Stern [3(3]. £rpa(z,0) = £ RP A(fc,0) is the static dielectric 
function. Wc have analytically evaluated the disorder-inclusive e(z,u) for which the results, presented below, appear 
to be new and we have utilized them in our numerical investigation. 

Let us recall the Lindhard (RPA) expression for the longitudinal dielectric function [lj|. In variables z and u and 
in 2D it reads as [301 



y 2 z r r 1 
£ RPA (z, u ,r) = i + / do 

TT JO JO 



qdq 



z 4 — (uz + iT + qz cos ( 



(2) 



= l + ^[F 1 (z,u,T)+iF 2 (z,u > T)} ) 



where we have split explicitly the DF £rpa {z, u, T) into the real and imaginary parts and have introduced the real 
functions (for real z and it) F\ (z, u, T) and F2 (z, u, T) as in the usual RPA expression of longitudinal dielectric function. 
Performing the q and 6 integrations in Eq. ([2]) [3l| we obtain, for a non-zero damping, 



Fi(z,u,r) = 2z+ - [Y- (z,u-) - Y- (z,u+)} 

z 

+ (u_ - l)Y + (z,u-) - (u+ - 1)1+ (z,u+) , 



(3) 



with u± = u ± z, 



F 2 (z,u,T) = -[Y + (z,U-)-Y + (z,u+)] 
z 

+ (u+ - 1) y_ (z, «+) - («_ - 1) y_ (2, «_) 



(4) 



Y ± (z,u) = ^= 



! z 2 ( u + iy + T 2 z 2 (u 2 -i) + r 2 



dz 



z 2 ( u -i) 2 + r 2 z 2 ( u -i) 2 + r 2 ' 



(5) 



In the case of vanishing damping (7—5-0 and r — > 0) the expressions ([T])-([5]) coincide with the Stern result [30l | with 

fi(z,u) = Fi(«,«,r)| r ^o = 2z + C-Ju 2 _ -\-C+Ju\-\, (6) 
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FIG. 1: (Color online) The contour plot of L(z,u) — Im[— l/e(z, u)] as a function of the variables z and u for r s = 2 and 
without (left panel) and with (right panel) damping (h^y = 15 eV). The dashed line in left panel shows the plasmon dispersion 
curve u r (z) (see the text for details) with 7 = 0. The numbers indicate the values of L(z,u). 



f 2 {z 7 u) = F 2 {z,u,T)\ r ^ a = D-^fl 



£Wl-<, 



(7) 



D ± = H (1 - \u ± \) , C±=H(\u±\-l) 



\u±\ 



(8) 



Here H(z) is the Heaviside unit-step function. The static DF involved in Eq. ([I]) can be found either from Eqs. ([3j> 
and (01 at the limits u — > 0, T — > or from Eqs. ©-© at u — > 0. The result reads 



£rpaM) = i+*-f(z) 



(9) 



with 



f (z) = -Mz,0) 



z, 
l 



z+Vz 



z > 1 



(10) 



To demonstrate the effect of the damping in Fig. [T] we show the contour plots of the energy loss function L(z, u) — 
Im[— l/e(z, u)] without (left panel) and with (right panel) damping. The plasmon dispersion function u r (z) in the 
left panel is also shown as a dashed line (the explicit derivation of the plasmon dispersion curve u r (z) without 
damping is given below in Sec. Ill Bl see Eqs. (fTI)]) . (fTT))) and Fig. [3]). The single-particle excitations energies huj sp = 
\hkvF ± fi 2 k 2 /2m\ (or u — \z±l\ in dimensionless units) are demonstrated as thick solid lines. As expected the 
energy loss function L(z, u) in the case of vanishing damping (left panel) is localized in the domains < u < 1 — z 
with < z < 1, \z — l|<u<z + las well as on the plasmon curve u r [z) where the function L(z,u) behaves as 
a Dirac (^-function and becomes infinity. In the case of non-zero damping (right panel) the energy loss function is 
broadened due to the damping and becomes non-zero also in the domains u < | z — 1 1 and u > z + 1 . 

Equations ((Tj)- (SJ) constitute the number-conserving DF for a 2D disordered electron gas. Deriving these expressions 
we have explicitly split the DF into real and imaginary parts assuming real variables z and u. An alternative (but 
equivalent) expression for this DF is derived in Appendix [A] which is valid for any complex u> and k. With this exact 
(within RPA and RTA) ex pres sion in Appendix [Al we then calculate the DF within small w, fc-approximation obtained 
previously in Refs. [271. l28l. 1291] and revisited in Appendix |A"1 The basic feature of this approximation is the prediction 
of the threshold condition for plasmon propagation which is absent in 3D (see, e. g ., R ef. J22l| ). Indeed the solution of 
the dispersion equation e (k,u) — 0, where e (k,u) is given by Eq. (|A8[) . reads [271 128|. |29| 



UJ r (k) 



1 + fcAxF 

1 + fcA TF /2 
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FIG. 2: (Color online) Real (Re[w r (fe)], left panel) and imaginary (Im[uj r (k)], right panel) parts (in a.u.) of exact (the lines 
without symbols) and approximate (the lines with symbols) solutions of dispersion equation vs k (in a.u.) for k ^ 2/cf, r s — 4, 
hj = 2.72 eV (solid lines), fvy = 8.16 eV (dashed lines), fi/y = 13.6 eV (dotted lines). 



where ui 2 (k) = 2itnQe 2 k/m is the plasma frequency for a 2D electron gas. The condition that uj r (k) has a real part 
(for plasmon propagation) leads to fc > fc*, where 



fc* — k> 



TP 



1 + 



7 



(12) 



with /ctf = 1/Atf = 2/arj. Thus, within small cj, fc-approximation, disorder in 2D electronic systems considerably 
softens plasmons; they cannot propagate for k < fc* and their dispersion relation is strongly altered relative to the 
collisionless case. However, since these results were obtained in small cu, fc-domain one can expect some modifications 
for large momentum transfers at k > fc*. Figure [2] shows the real (left panel) and imaginary (right panel) parts of the 
solutions of the dispersion equations with approximate (Eq. (p~Tj) ) and exact dielectric functions, Eqs. (| A3|l - (|A6[) . For 
simplicity we consider the case fc ^ 2fci? when the function Q in Eq. (|A5p vanishes. Note that the condition fc > fc* 
together with the inequality fc < 2fci? requires that H^/Eyi < 4rJ 2 v / l-414r s + 1 with En = me 4 /h 2 ~ 27.2 eV. It is 
seen that the slope of the imaginary part of u r (fc) (right panel) is dramatically changed at some value of fc where 
the expression under square root in Eq. (jTTJ) changes the sign. For small u>, fc-approximation this value of fc is given 
by Eq. (TT2"]) . As pointed out in Appendix [5] the approximation (|lip is valid when one neglects the single-particle 
energy hwk = h 2 k 2 /2m with respect to HkvF- Therefore, in general, we expect good agreement between approximate 
and exact to r (fc) for small momentum fc, as shown in Fig. [5] However, with increasing 7 the approximate dispersion 
relation (fiT|) fails to predict io r (fc) correctly. As shown in Fig. [2] (left panel, dotted curve with hj — 13.6 eV) at small 
fc (fc < fc*) the energy of plasmons uj t (fc) is not exactly zero as predicted by Eq. (fil"]) although the probability of 
plasmon generation is strongly reduced due to the relation Re(av) -C Im(oj r ). Moreover, in contrast to the predictions 
of approximation (|11[) in this case with increasing momentum fc the real part of uj r vanishes and plasmons cannot 
propagate any more. 



B. Strongly coupled electron gas: Beyond RPA 

In this section we consider exchange-correlation interaction effects via local-field corrected (LFC) DF but we neglect 
the disorder (i.e. 7 = 0). To include disorder in a fully interacting electron gas at a microscopic level is rather involved, 
and no analytical calculations of e(fc, u>) without restrictions on fc and lo are still available. An attempt to involve strong 
correlations in RTA and within small fc, w-approximation (see Eq. pT|) ) has been done in Ref. [28[. Instead we employ 
here the LFC dielectric function and demonstrate some useful results which have not been considered previously. Our 
discussion below is based on the LFC dielectric function of a fully DEG see, e.g., Ref. |§| (in dimensionless variables 
z and u) 

e(z,u) = l+ ffi") (13) 
1 — G(z)T'(z 7 u) 
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TABLE I: The critical wave numbers (dimensionless) z c and the minimum values A c = «i-(2min) of the dispersion function u r (z) 
for some values of the density parameter r a . z c °^ and A^ ' represent the smae quantities but for non-interacting 2D electron 
gas. 
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Zc 


0.126 


0.288 


0.390 


0.457 


0.507 


0.546 


0.579 
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0.629 


zi 0) 


0.135 


0.345 


0.510 


0.638 


0.748 


0.845 


0.932 


1.013 


1.089 


\c 


1.116 


1.264 


1.358 


1.421 


1.469 


1.507 


1.538 


1.565 


1.588 


xi 0) 


1.122 


1.304 


1.440 


1.543 


1.629 


1.704 


1.770 


1.831 


1.886 



where V(z, u) is the polarizability of the free-electron gas obtained in RPA by Stern [3(| 

x 2 

V{z,u) = erpa 0,«) - 1 = ih( z > u ) + ih{z,u)\ (14) 

with £rpa(z,m) = £rpa(z,w, T ~ > 0), where £rpa (z, u,T), /i(z,it) and f2(z,u) are given by Eqs. ([2]), (O and__ 
respectively. Note that our definition of the functions fi(z, u) and f2{z, u) differs from the definition given in Refs. 
by a factor of —1/2. G(z) is the LFC function, which includes the effects of exchange-correlation interactions. Within 
a sum-rule version of the self-consistent approach, Gold and Calmels presented [32[ a parameterized expression G(z) 
for the 2D electron gas, 

G{z) = -= ( ZG ^l 2 , ■ (15) 
^G{ 2 (r s ) + Z 2Gj 2 (r s ) 

1 /3 

The coefficients Go (r s ), G\2{r s ) andG22(?~s) are determined by Go (r s ) — 1.983r s ' , Gi2(r s ) = 1.626Gi2(r s ), G22(r s ) = 
\[2r a 1 ^C22{r s ), with Gi2(r s ) = air] 1 , C22(r s ) = c^g 2 , and the parameters a%, ct2 and 71, 72 can be found in 
Ref. H. 

Now we consider the exact solution of the dispersion equation e{z, u) — for an interacting electron gas when the 
DF is given by LFC expression (|13p . From Eqs. (O-© and (fH?)) . (jT^J) it is seen that the collective plasma modes 
(plasmons) can propagate with the frequency and momentum u> and k (or u and z) which lie in the domain u ^ z + 1 
where f2(z, u) = and Im[e(z, u)] — 0. In this domain the dispersion equation has an exact analytical solution which, 
in Lindhard's dimensionless variables, is given by 

u 2 r (z) = l + z 2 [azg(z) + l} 2 + 1 (16) 

azg (z) [azg (z) + 2\ 

with a = V2/r s and g (z) = [1 — G(z)]~ 1 . It is straightforward to check that the solution (fL6|) indeed satisfies the 
condition u r (z) ^ z + 1 for arbitrary z. However, an inspection of the dispersion equation shows that this solution 
exists only for the wave numbers from the domain ^ k ^ k c (or ^ z ^ z c ) where the critical wave number z c 
is obtained from an equation u r (z c ) = 1 + z c , i.e. in this point the plasmon curve u r (z) touches to the boundary of 
the single-particle continuum u = 1 + z. Explicitly, the critical wave numbers are determined from transcendental 
equation 

az\g (z c ) [2 + az c g (z c )] = 1. (17) 

Table |T] shows the quantity z c and the minimum of the dispersion function A c = u r (z ln - m ) with u' r (z m i n ) = for 

some values of the density parameter r s . The critical wave numbers and the quantities A c (labeled as zi ^ and A^ ' , 
respectively) are also shown for non-interacting electron gas, i.e. with G(z) — and g(z) = 1. These quantities are 
important for evaluation of the SP in Sec. IIII1 

We can present the dispersion expression (|16[) obtained above, in the usual form 



^ {k) = ^ {k) 2 + kX TF g(k) (18) 

2 2-l g (fc)+fcA TFg (fc) ^ ft 2 fc 4 a 
+k Vf 2 + k\ TF g(k) + 4^ [1 + " Atf5 (fc)] 
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Z, 

FIG. 3: (Color online) The dispersion curve u r (z) for free (lines with symbols) and interacting (lines without symbols) electron 
gas vs z for r s — 0.5 (dashed lines) and r B — 4.0 (solid lines). The oblique solid line corresponds to uj = kvF + hk 2 /2m (or 
u = z + 1). 

which for vanishing exchange-correlation interactions (i.e. at g (k) = 1) reads 

.1 (*) = .1 (k) + *V,f±|j^ + % (1 + fcA TF) 2 . (19) 

This exact (within the employed model) dispersion relation may be compared with an approximate result derived by 
Fetter within a hydrodynamical approach [331 ] - Equation (I19p agrees with the hydrodynamic result if the last term 
(the single-particle energy) in this expression is neglected and the coefficient at k 2 v 2 F is replaced by a constant factor 
1/2. It should be emphasized that in general and at long wavelengths ui r (k) from Eq. (JTSJ) for an interacting 2D 
electron system varies like k 1 / 2 independently of the LFC G(fc) and in contrast to the 3D case. This latter behavior 
seems first to have been suggested by Ferrell [lij and later investigated in more detail by Stern [3(| (see also the 
review paper [35[). It arises from the electromagnetic fields in the vacuum surrounding the plane, with an associated 
reduction in the screening. Since uu r (k) increases monotonically from zero, an external perturbation of arbitrarily low 
frequency can always excite collective modes. Hence, the characteristic 3D absorption edge at constant 3D w p is here 
entirely absent. Moreover, the group and phase velocities both diverge like k^ 1 / 2 as k — > 0. 

Figure [3] shows the plasmon dispersion curve u r (z) for interacting (the lines without symbols) and non- interacting 
(the lines with symbols) electron gas, i.e. Eqs. HU) and p^|) . respectively. The points where the plasmon curves touch 
the single-particle excitation boundary are given by z c or , see Table |U It is seen that the exchange-correlation 
interaction may strongly reduce the values of u r {z). It must be pointed out a technical but important detail which, 
to our knowledge, has not been yet discussed in the literature. From Fig. [3] it is seen that in u, z plane the plasmon 
curve u r (z) has a minimum which is absent in usual units u,k where ui r {k) is a monotonic increasing function. By 
interchanging the z and u axes in Fig. [3] one obtains the plasmon dispersion curve z r (u) which, however, in contrast 
to the 3D case has two different branches with increasing (z r i(u)) and decreasing (z r 2(u)) dispersion functions (at 
the minimum of u r (z) both z r i(u) and z r 2(u) curves contact each other). The dispersion relations z r i(u) and z r 2(u) 
can be provisionally treated as the "single-particle" and "plasmonic" relations, respectively. Therefore when one 
attempts to perform ^-integration in Eq. (|20p before u- integration, as was done in Ref. [||, the double integration in 
the SP is reduced to two line integrations along the contours z r i(u) and z r 2(u) and both of them contribute to the 
SP. In other words in this case the energy loss function L(z,u) introduced above contains two Dirac ^-functions. In 
fact, we see from our numerical calculations that near the SP maximum the contribution of z r \{u) is not necessarily 
small compared to the contribution of the other one, z r 2{u). Although the total contributions of both in the SP are 
in general much smaller than the purely single-particle contributions. This is a violation of the Lindhard-Winther 
equipartition sum rule [17| which we further discuss in Sec. IIII Al To avoid this technical problem in the numerical 
calculations it is easier to perform first in Eq. (|20p the it- integration and then using the dispersion function u r {z) 
given by Eq. (frp|) . 
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III. STOPPING POWER 



With the theoretical formalism presented so far, we now take up the main topic of this paper. This is to study 
the stopping power (SP) of a point-like ion in a 2D degenerate electron gas as well as to show how collective and 
single-particle excitations in the target medium DEG contribute to the SP. And, as in the previous section, we shall 
present new theoretical results within the linear response approach. We consider two models for a DEG in 2D. 
(i) A disordered DEG for which we use a number-conserving DF given in Eqs. (HJ)-(j5l) . For this case we present 
analytical calculations and new results for the SP in a low- velocity limit, (ii) A strongly coupled DEG with a DF 
which includes LFC, Eqs. (fT5| and (|14[) . This case has been studied in Refs. [f| 0] where the leading term in a 
high-velocity limit of the SP is calculated using a plasmon-pole approximation. This calculation is supported by a 
more rigorous treatment, again for the leading term only, in Ref. |8] which is based on a method of moments and 
includes electron-electron interactions. Now the leading term happens not to depend on electron-electron interaction. 
It is then of interest to calculate analytically the next non-vanishing terms of the high-velocity SP. As shown below 
these terms are significantly modified by electron-electron interaction and thus are more involved than the leading 
term. 

We consider an external point-like projectile of carge Ze moving with velocity v in a homogeneous and isotropic 
2D electron medium characterized by the dielectric function e(fc, lu) or s{z^u). Then in the linear response theory the 
SP which is the energy loss per unit length by this projectile is given by [3J, l|| 



8£ Z 2 r udu 



Im— rzdz, (20) 



7i"X 4 A Jo V\ 2 - u 2 Jo e{z,u) 

Here A = v/vp, So = e 2 /a§ ~ 5.132 GeV/cm = 51.32 eV/A. We have used the Lindhard variables z and u introduced 
in Sec. HH In our calculations we shall consider the range of v for which the linear response theory is found to be 
adequate [361 ]. 



A. Equipartition sum rule 



With the theoretical formalism presented so far, we now take up one of the main topic of this paper. This is to 
study how collective and single-particle excitations in the 2D electron gas contribute to the SP. This problem was 
first addressed by Lindhard and Winther [ijj (LW) for a 3D degenerate electron gas without damping (7 = 0). They 
formulated an equipartition sum rule which states that an integral proportional to that in Eq. (|20[) 



3 (u) — 3 sp (u) + 3 p (u) = I Im— -zdz (21) 







receives equal contributions from plasmon (Qp) (with 0<z<u— 1) and from single-particle excitations (3 sp ) (with 
u — 1 < z < u + 1), respectively. The functions 5 p (u) and ^s sp (u) may then be written as 



% («) = ^ Im-7- —,zdz = ^ y , (22) 
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-is 



(«) 



Im 



e (z, u) 



zdz. 



(23) 



Here z r (u) is the solution of the dispersion equation e{z, u) — (the inverse of the dispersion function u r (z)). This 
equipartition rule is valid for sufficiently large u, u > Um, where the threshold value u m in 3D case is obtained from 
the equation z r (u m ) — u m — 1. In recent works (20l . [23l.l37| we have shown that the LW equipartition rule does not 
necessarily hold for an extended charged projectile e.g. a diproton cluster in a 3D degenerate electron gas without 
disorder (7 = 0) as well as for a point-like ion in a disordered DEG. We have established some generalized stopping 
power sum rules. In this section we briefly show that the LW equipartition rule is also violated for a 2D electron gas. 
In the present context it should be emphasized that the plasmon contribution given by Eq. (|22[) contains indeed two 
terms, with z r \ (u) and z r 2 (u), as discussed above. The existence of both branches requires the threshold condition 
u > u m , where u m is the minimum value of the dispersion function u r (z) shown, e.g. in Fig. [3J However, it is 
clear that the contribution of z r \ (it) vanishes at u > u ci where u c — u r (z c ) > u m (the point where the plasmon 
curve touches to the single-particle excitations boundary). For simplicity we consider below only the domain u > u c 
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where only z r 2 (it) = z r (u) contributes to the SP integral (j22)) . As an example we employ the DF (fT3)) together with 
Eqs. ©-(jSJ) and (jT5J) for an interacting DEG. The simplest way to show the violation of the LW equipartition rule 
in 2D is to calculate the asymptotic values of the contributions 3 p (u) and 3 sp (u) at u > 1. The inverse dispersion 
function z r (it) for 2D interacting DEG is evaluated in Appendix [Bl see Eqs. (|B5I) - (|B6I) . Using these expressions it is 
straightforward to calculate the single-particle and collective contributions to the SP integral which at u ^S> 1 become 



4a 2 u 4 1 2u 2 4u 4 



29 

h 2 (r s ) -3*r(r s ) + — 



sp ^ ^ Aau i 4u 2 2au 3 



1 



(r«) 



1 

8m 1 



(24) 
(25) 



Here xq (r s ) and x(r s ) are defined in Appendix [Bl From the above expressions it is clear that the contribution of 
the collective excitations is much smaller than the contribution from single-particle excitations, 3 P (u) <C 3 S p (u 
which indicates the violation of the LW equipartition rule. A similar result has been found numerically in Ref. 
and is supported by our own numerical calculations. Of course, Eqs. ([M]) and (f2"5]l are not strong results. An exact 
treatment can be developed on the basis of the integration contour on the complex z-plane suggested by LW [l7j and 
investigated in details in Ref. [23[ . The technique developed in [23[ is independent of the dimensionality of electron 
gas but requires a necessary analytic continuation of the DF in the complex z-plane, that is e(— z*,u) — e*(z,u), 
where the asterix indicates a complex conjugate quantity. It is easy to see that this condition is violated for a 2D 
electron gas. For simplicity let us consider non-interacting DEG with the DF given by Eq. §2§ in the integral form 
and with T — > +0. In this case one can easily check that £rpa (— z* , u) — 2 — £ RP a (z, u) (a similar equation can be 
obtained for an interacting electron gas) . Therefore an arbitrary function of the form 

C 

e cS (z, u) = 1 + - [erpa (z, u) - 1] (26) 

with an arbitrary constant C defines an effective DF of a 2D electron gas which satisfies the required condition, i.e. 
£cS (— z* , u) = e* ff (z, u). Applying now the contour integration technique developed in Ref. 23] one can strongly prove 
that the single-particle and collective excitations contribute equally to the SP integral (l2Tj) where the DF e(z,u) is 
replaced by the effective one, £ e g(z, u), given by Eq. (|26|) . Thus the LW equipartition rule holds also in 2D treating 
the effective DF instead of e(z, u). In this case it is straightforward to check that at u 3> 1 the leading order terms 
of the collective and single-particle excitations are given by 3 P (u) — 3 sp (u) ~ 7r/(4cm 2 ). The physical origin of the 
modification of the equipartition rule in 2D is the change of the nature of the Coulomb potential (in Fourier space it 
behaves as ~ 1/fc in 2D) and as a consequence the long-wavelength dispersion relation: the plasma frequency behaves 
as ~ A; 1 / 2 in this limit. Technically this modification introduces an extra non-compensated z variable as a prefactor 
in Eq. first line, which changes the analytical properties of the DF. Introducing an effective DF (1^)) we formaly 
replace the 2D Coulomb potential by the 3D one without affecting the polarizability of the 2D system. This recovers 
formally the 3D-type dispersion relation with constant plasma frequency and hence the equpartition rule. 

B. Low- velocity limit 

Let us consider SP for slow projectiles, with v -C Dp. A consequence of the 3D linear response theory, confirmed 
by experiments, is that for ion velocities v low compared to the Fermi velocity vp, the stopping power is proportional 
to v (see, e.g., the latest experiment |38j). The coefficient of proportionality may be called a friction coefficient. A 
similar linear behavior of the SP, S ~ v, is expected in 2D case Refs. 0, 0, [E S 13 • Using analytical results obtained 
for £rpa(z, u, r) the general expressions for SP follow from Eqs. (I20p -((5l): 

X Jo [z 2 + xfl z )\ X 

where the dimensionless friction coefficient 5i(r,x 2 ) depends on the target properties and hence also on the dimen- 
sionless damping parameter T. We have introduced the following functions 

i / M p/M-^,,r)l t (28) 
r ip{z,r) 

1>{z, T) = F 1 (z, 0, T) = 2z + - z [*_ (z) - (-z)] -(z + 1) $+ (z) - (z - 1) $ + (-z) , (29) 
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(z-l) 2 +r 2 z 2 (z 2 -l)+T 2 



y/2 \| V z 2 (z + 1) 2 + r 2 ± z 2 (z + i f + r 2 



(30) 



The static screening function / (z) is determined from Eq. (fit) 
When the damping vanishes (r — > 0) Eq. (|29j) becomes 



V(z,r)->2/(z) 



2r 



:ff(i-z) + o(r 2 ) 



where (z) is the Heaviside unit-step function. Therefore 



E(z,T) 



and from Eq. (|2Tjl we find 



r-»o 



=H(l-z) 



(31) 



(32) 



®( r .X 2 )lr^o 



% (x 2 ) 



, X 



2 1-x 4 



z dz 



(33) 



2(X 4 - 2) 



S(x 2 ) 



with 



g(x) = 




X > 1 



a; < 1 



(34) 



The last expressions (|33|) and ([34]) are known results derived previously within RPA in Refs. 0, @| • Interestingly, in 
a low- velocity limit this SP completely agrees with the result obtained within a binary collision approach Ref . 9] . In 
left panel of Fig. Q] we show the ratio of the disorder-inclusive friction coefficient 5R(r,x 2 ) and 5Ro (x 2 ) vs damping 
parameter fi/y for two values of the density parameter r s = 1 and r s = 2. To gain more insight in right panel of Fig. [4] 
we show the friction coefficient 5R(r, x 2 ) vs r s for some values of the damping parameter 7. As expected, the friction 
coefficient and hence the SP at low velocities increase with an increasing damping parameter 7; this was previously 
reported for 3D in Refs. [13, H^HJ]. The behavior of 5R(r,% 2 ) at fixed 7 and at increasing density parameter r s is 
particularly noteworthy. At small damping the friction coefficient decays monotonically with r s while at large 7 it 
may also increase for large r s . We will further discuss this behavior in Sec. IIVI 

The approximation (|27p implies that the SP is proportional to velocity. The velocity region in which the linear 
proportionality between SP and the projectile velocity holds may be inferred from the numerical calculations (see 
Sec. IIV[) . It is seen from those results that the approximation (|2"T]) remains quite accurate even when A becomes as 
large as ~ 1. 



C. High- velocity limit 

Consider next the limit of large projectile velocities in the case of strongly interacting DEG with the dielectric 
function Eqs. (fT51) - ([T51) . In this limit the general expression (f2"0"l) forpoint-like projectiles with charge Z moving in 
either interacting or free electron gas reduces to the simple formula [3| 

ttSqZ 2 27r 2 n Z 2 e 4 

S ~ — — — = (35) 

X A nv 

which does not contain the gas electron mass m anymore; m and also the effects of electron-electron interactions 
appear only in the higher terms of the expansion. The other main discrepancy between the 2D and the 3D results 
is that the stopping power decreases as 1/v instead of behaving as \n(v)/v 2 in the 3D case. In the presence of 
interactions the next order terms are shown to be significantly modified. We derive below a generalized expression 
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FIG. 4: Left panel, the ratio of the friction coefficients with and without damping r](T,x 2 ) = X 2 )/%(x 2 ) vs damping 
parameter tvy (in eV) for various densities. r s = 1 (solid line), r a = 2 (dashed line). Right panel, the friction coefficient K(r, x 2 ) 
vs density parameter r s for 7 = (solid line), hf = 1 eV (dashed line), Try = 10 eV (dotted line) and ftrj — 20 eV (dash-dotted 
line) . 



for SP, in a high-velocity limit, for point-like ions. In order to show how SP in a high-velocity limit is affected we 
consider expression (|20|) rewritten as follows: 



S = 



2Z Z 2 f x A(u)du 



— ' 



where 



A(u) 



4a 



u5 (u) 



4a 



Im 



e (z, u) 



zdz 



(36) 



(37) 



and 3 (it) is the total contribution of the collective and single-particle excitations to the SP integral defined in Sec. lIII Al 
(see Eqs. ([2"Tj) - (j2"3")) '). For further progress it is imperative to calculate the asymptotic behavior of the function A (it) 
at u — > co. For collective and single-particle excitations these asymptotic forms are given by Eqs. (I24|) and l|25[). 
respectively. Using these expressions we arrive at 



A(«) 



1 



C2 C 3 , _ 4 s 

-T + — r + O [U 



(38) 



for u — * oo and with the expansion coefficients 

C 2 = 



c. 



5- [1 - xt (r,)] 



(39) 



Here the parameter (r s ) depends on the exchange-correlation interactions and is given explicitly in Appendix IB] 

Below we calculate the SP up to the order 0(v~ 4 ) thus neglecting the terms with O (A -5 ). First the SP (l36|) can 
be represented in the equivalent form 



S = 



ttS Z 2 
ttA 3 



hi 

1 + T 
i 

+ 4^ 



2ttA 2 



1 - 



1 



A 



A + VA 2- 1 

which is convenient for further calculations. Here hi is a constant 

2 



A + VA 2 - 1 

In ! A + ^A 2 



- $ 2 (A) 



(40) 



$i (A) 



h = 



A (w) - 1] du 



(41) 
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and the other quantities are function of the ion velocity: 

- 1 ) [A (u) - 1] du 

A(u) 



~u 2 /\ 2 

1 



2A 2 r A / 

'CWi I y/l-v?/\ 2 



1 



1 



9i 



C 3 



du 



(42) 



1 



$ 2 (A) = 4A J [A («) - 1] du. 



(43) 



For the derivation of Eq. (|4"0"1) we have used some elementary integrals [3l|. In Appendix [B] we prove that hi = 0, see 
Eq. (|B8|) . This relation can be regarded as another SP sum rule for an interacting DEG. 

For a calculation of the SP up to fourth order v~ 4 we need the asymptotic behavior of $2 (A) up to the first order 
(v^ 1 ) which can be obtained from Eqs. (|3"8|) and (|43[) 

(44) 



$ 2 (A) = l + ^ + 0(A- 2 ), 
and only the leading term of $1 (A). We denote this leading term by $1 (A)| ; 



In /12 and using Eq. (I42|) we obtain 



In /19 = — — 

C 3 



[A (u) - 1] u 2 du + 



A (u) - 1 - 

qiA I/O 



2 du- - \ - -. 



(45) 



The coefficient In /12 is explicitly evaluated in Appendix [B] and entirely depends on the density parameter r s , see 
Eq. (|B14p . Thus substituting Eqs. (JSJ) and (gSJ) into (gDJ) and setting hi = we finally obtain 



S 



irZ Z 2 



ttA 3 



ln(2/i 2 A) 



(46) 



It is seen that in the correction term (the second term in Eq. (|46|) ) the mass of electron enters through the Fermi 
velocity vf — hkp/m. A limit to the non-interacting DEG is performed by taking the limit h\ (r B ) — > 0, i.e. setting 
C3 = 3/2a = (3\/2/4)r s (see Eq. (f39| ). In this limit the coefficient ln/12 is given by Eq. (|B15|) . In the general case 
of non-vanishing exchange-correlation interactions it is too difficult to draw some conclusions from Eq. (14"6"|) about 
how these interactions affect the high-velocity SP. Numerical calculations of Refs. @, 0] show that these interactions 
strongly increase the SP up to the intermediate velocity range with v ~ vf- We support this conclusion by our own 
calculations (not shown here) which also indicate that the asymptotic SP (|4"6"|) remains quite accurate also in the 
intermediate velocity range. 

We close this section with the following two remarks. First, the high-velocity SP Eq. (|4"6")) is also valid for a general 
LFC function G(k). The derivations above and in Appendix IB1 show that only the asymptotic values of G(k) at 
k — > 00 and k — > contribute to Eq. (|4*6^1 . At short wavelengths G(k — > 00) = Goc^s) is constant (see, e.g., Ref. |32|). 
At long wavelengths the LFC function behaves as G(k — > 0) ~ K(r s )k/kp, where the constant fc(r s ) is related to the 
compressibility of a 2D electron gas through compressibility sum rule. The latter for a 3D electron gas is discussed in 
[39l | . and for a 2D electron gas in [4(j. Thus in the general case of arbitrary G(k) the quantities xi(r s ) and X2(r s ) in 
Eq. (|4"6")l are replaced by xi(r s ) = {\ / i)G x {r s ) and X2(r s ) = G ^(r s ) / 'Ak(t s ) , respectively. Second, a similar procedure 
is applicable to evaluate the high-velocity corrections also for a disordered 2D electron gas. While the high-velocity 
SP (|46|) does not contain the terms of the second v~ 2 and third v~ 3 orders, some preliminary investigations by us 
show that for a disordered DEG this SP involves also the terms of the order B\v~ 2 , B2V~ 2 hxv and B%v~ 2 In 2 v, where 
the constants B\, B2 and B3 depend on 7. Therefore the corrections to the high- velocity SP would be much more 
sensitive to the ion velocity than those predicted by Eq. (|4*6")l . 



IV. NUMERICAL CALCULATIONS 



Using the theoretical results obtained in Sees. HT1 and IIII1 we present here the results of our numerical calculations of 
stopping power for a 2D target material with the wide range of the density parameter, 0.1 ^ r s ^ 5. The parameter r s 
varies from the small (free DEG) up to the large (strongly interacting DEG) values. As examples of 2D target material 
we have considered two models. An interacting DEG whose linear response function includes the exchange-correlation 
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FIG. 5: Left panel, the SP (in a.u. and in logarithmical scale) of a proton vs v/vf moving in an interacting electron gas for 
r s —0.1 (dotted line), r s = 2 (dashed line), r a — 5 (solid line). The lines with and without symbols correspond to the numerical 
evaluation of Eq. (|20p with Eqs. (|13p - (|15[) and asymptotic expression (|46|l . respectively. Right panel, the dashed and solid lines 
(with symbols) from the left panel in linear scale. 



effects via LFC and is given by Eqs. (fT3")) - ([T5|) . This case has been investigated previously in Refs. @, 0|- In Fig. O 
left panel we compare the exact (the lines with symbols) and asymptotic (the lines without symbols) SPs calculated 
from Eqs. ([20| . (|13jl - (fT5|) and (|46|) . respectively. It is seen that the asymptotic expression (j46| is very accurate and 
at v > vf practically coincides with exact SP. In general we have found that the higher order correction in Eq. (J46J) 
(the second term) is small compared to the leading term. However, the role of this term becomes more and more 
pronounced with increasing the density parameter r s , i.e. with increasing the exchange-correlation interactions. We 
have also compared our numerical calculations with the results obtained by Wang and Ma Two major differences 

have been found. First, the LFC dielectric function (fTB"]) for a fully degenerate electron gas predicts a threshold ion 
velocity for plasmon excitations. In view of the discussion in Sec. Ill Bl the plasmons are excited at A > A c , where the 
critical (dimensionless) velocity A c is the minimum value of the dispersion function u r (z), Eq. (|16[) . and can be found 
from the equations A c = u r (z m i n ) with u' r (z m i n ) — (see also Table [TJ) . The velocity threshold changes sufficiently 
the slope of the SP and at A = A c one expects a characteristic discontinuity of the derivative of the SP (the SP itself 
remains naturally continous at A = A c ). In contrast to Refs. [1,0] this feature is clearly visible in Fig. [51 right panel 
(see also the solid lines in Figs. [6] and [7]). Such behavior of the SP at A = A c has been observed previously in 3D (see, 
e.g., an d references therein). 

Second, we have found that for the same conditions (i.e. for the same r s ) the SP in our case is considerably smaller 
near maximum than those obtained in Ref. 6;]. Moreover, there is no agreement between the results obtained in 
Refs. @ and Q, e.g. for r s = 1 and r s = 5, where Ref. [(| predicts in whole velocity range much larger SP than the 
latter. Apparently this is because the polarizability of the free-electron gas employed in Ref. @ somewhat differs from 
original expression derived by Stern [301 ( see a l so Eq. (fTij) with Eqs. ©, GO; the algebraic square roots in Eqs. © 
and ||7J) are missing in Ref. [fjj. These square roots are recovered in Ref. 7| but nevertheless one of two plasmon 
branches is ignored as discussed in Sec. Ill B I which may yield smaller value of the SP. 

Within the second model the target material is modelled as an electron gas whose linear response function, within 
RTA, is given by Eqs. (!])-© with 7 as a model damping parameter. In Figs. [5] and [7] based on numerical calculations 
we choose four values of 7: fry = (solid lines), fry = 1 eV (dashed lines), fry = 10 eV (dash-dotted lines), and 
fry = 20 eV (dotted lines). The values < fry < 1 eV are comparable with the damping parameters (related inversely 
to the collision times) in some 3D metal targets, e.g., Al for which fry can be ~ 0.1 eV. The last values fry = 10 and 20 
eV correspond to the damping parameter, e.g. in carbon. It is seen from Figs. [5] and [7] that the SP is broadened with 
increasing damping (i.e. with increasing 7) and this effects is more pronounced for small densities (i.e. at large r s ). 
Of course, the value r s = 4 in Fig. [7J is somewhat far beyond the RPA employed for deriving the dielectric function 
Eqs. (d])-©. However, treating this case as a qualitative example we look for some complementary information about 
the effect of disorder at large r s . In particular, at either vanishing or small damping with 7 = and fry = 1 eV, 
respectively, Fig. [JJ predicts a modification of the linear friction law (see Eq. (|27[) ) which now approximately behaves 
as ~ v 3 . This is v 3 law obtained e.g. in Ref. [13] within linear response theory for classical 3D plasma and supported 
by the numerical simulations [36[ . According to Fig. 0] (right panel) at small 7 the linear friction coefficient decays 
with r s and may be smaller than the cubic friction coefficient (~ v 3 ). However at strong damping the linear friction 
coefficient again becomes dominant and the SP at small velocities behaves as ~ v (cf. Figs. 2] (right panel) and[7J. 
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FIG. 7: Same as in Fig. [6] but for r a = 4. 



V. SUMMARY AND CONCLUDING REMARKS 



In this paper we have presented a theoretical study of the stopping power of point ion projectile in a degenerate 2D 
electron gas. The later is modelled within two different approaches namely (i) as a system containing disorder due to, 
e.g., electron-impurity interactions and (ii) including exchange-correlation interactions of the electrons. In the course 
of this study we have also derived some analytical results for the disorder-inclusive RPA linear response function and 
for the corresponding plasmon dispersion relations. These analytical results go beyond those obtained previously in 
Refs. (27l.[28l.[29| within small-A:, u) approximation. Also for the model (ii) we have found an exact dispersion relations. 
After a general introduction to the SP of an ion in Sec. [TTI theoretical calculations of SP based on the linear response 
theory and using the models (i) and (ii) are discussed in Sec. Mil A number of limiting and asymptotic regimes of low- 
and high-velocities and vanishing damping have been studied. These approximate expressions are well supported by 
our numerical calculations. Special attention has been paid to the equipartition sum rule in 2D. In Sec. Mil employing 
the model (ii), i.e. LFC dielectric function for an interacting DEG, we have shown that this rule does not necessarily 
hold in 2D and may be satisfied introducing an effective dielectric function (|26|) . The theoretical expressions for a 
number of physical quantities derived in this paper lead to a detailed presentation, in Sees. ITIllIVI of a collection of 
data through figures on SP, friction coefficient and the dispersion relations. For the damping parameter, we have 
chosen a wide range of values ^ hrf 20 eV; the damping parameters for some 3D metal and semiconductor targets 
fall within this range. The results we have presented demonstrate that with regard to several physical quantities of 
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primary interest the difference between RTA and usual RPA without damping is significant. 

It is of particular interest to study the high- velocity limit for the SP of an ion beam. Such asymptotic expressions 
contain some useful information on a projectile ion structure factor and specially on the target medium properties. 
Eq. (|4"rJ)) with Eq. (|B14p which are a generalization of the asymptotic formula obtained in Refs. [1, 0, [H, Q can be 
used for analyses of experimental data on high-energy beam interactions with 2D target material. We note that the 
analytical method developed here for the derivation of high- velocity SP is general and may be applied within a linear 
response treatment for other types of projectiles, e.g. extended multicharged ions, as well as for any particular form 
of the linear response function e (z, u) for the target material. For given target material this approach requires only 
the asymptotic form of the plasmon dispersion relation at high u = uj/kvF and the frequency moments of the energy 
loss function. For a disordered DEG (model (i)), however, some modifications occur when one includes the damping 
in the DF. For instance, at large frequencies the energy loss function Im[— l/e(z, u)] for a disordered DEG behaves 
as joj 2 (k)/oj 3 and obviously the third frequency moment of this function does not exist, see Eq. (|B9[) . This requires 
some additional investigation of the third moment sum rule for this case which in turn is important for evaluation of 
high-velocity SP, see Sec. [TTTCl 

We shall make some brief remarks on the RTA in the linear response function. In the present study the disorder- 
inclusive linear response function containing in the RTA has been considered only in RPA. Going beyond RPA with 
electron-electron interaction and disorder treated at the same microscopic level is a difficult task. We may mention 
that recently the linear response function in 3D has been considered in RTA which conserves the particle number, 
momentum and energy [43l . |44| (see also references therein) . We intend to extend this model with fully conserving 
(number, momentum and energy) linear response function for 2D electron gas. 

In our calculations of SP and related quantities we have modelled the disordered 2D target medium as an electron 
gas whose linear response function is constructed in RTA in order to include scattering of electrons with disorder 
impurities. The numerical values of the phenomenological quantity 7 used in our calculations are within a physically 
expected range for the specific target medium. In principle 7 can be calculated to varying degrees of approximations. 
In the simplest approximation, its inverse can be calculated through Fermi's golden rule for a model electron-impurity 
potential. This may allow us to see how SP and related quantities depend on the target properties through their 
influence on 7. 

We expect our theoretical findings to be useful in experimental investigations of ion beam energy losses in solids. 
One of the improvements of our model will be to include some short-range correlation in the linear response function. 
Another interesting issue not considered here in detaile is the effective DF (|2l))) for 2D interacting DEG. Our goal is to 
find physical motivation and basis for this type of DF. A study of this and other aspects will be reported elsewhere. 
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APPENDIX A: DIELECTRIC FUNCTION OF DISORDERED ELECTRON GAS 



In this Appendix we give an alternative derivation of the disorder-inclusive DF which is valid in the entire complex 
u>, fc-plane. Since we are going to compare our results with previous derivations in Refs. (27l . [28l . [29| here we use the 
usual energy (o>) and momentum (k) variables. Performing the q and 9 integrations in Eq. ^ without splitting this 
expression into real and imaginary parts, for arbitrary u> and k complex variables we obtain 



£rpa (k,uj) = 1 + - — 

KQq 



2lu, 



2 vl + 



\J (uj + + uj k y 



k 2 Vp 



(Al) 



where ui + = uj + 17, u)k = hk 2 /2m. Here huk is the single-particle energy and Atf = Qo/2 plays a role of the Thomas- 
Fermi screening length which is constant in 2D case [33| . The multi- valued functions in Eq. (|A1|) must be understood 
in the following way. (i) The imaginary parts of the square roots are positive, (ii) The signs of the real parts of the 
square roots with u + ± ujk are taken with the sign of the expressions 



|u±| \uj±hk 2 /2m\ 



u± lu ± hk 2 /2m 



(A2) 
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These two conditions completely fix uniquely the values of the square roots entered in Eq. 
conserving DF is now evaluated using Mermin-Das formula, Eq. JT]) 



Here 



1 - 



w + ijQ (k, uS) 



ka [' P (k, cj) + ry [Q (k, u) - 1] 

\J (UJ + - U! k ) 2 - k 2 V 2 p + \J (UJ + + LO k f - k 2 v 2 F 



Tl)- The full number- 



(A3) 



(A4) 



Q(k,u) = -^- — [P{k,w)-wi 



v{k) 



£RPA(fc,0) 
£RPA(fc,0) - 1 



= 1 



v(k)-l 



kao 



kao k/2kp 
~f(k/2k F ) 



(A5) 



(A6) 



and the function / (z) has been introduced by Eq. (jTOJ) . 

Now let us consider the limit of small momentum-energy transfers, i.e. we assume that k -C 2k p and hu> -C Ep- In 
this case /(z) = z and the function Q in Eq. (|A5[) vanishes. In addition neglecting the single-particle energies hw k in 
Eqs. (lATj) . R3|l and (O we obtain 



£RPA (fc, w )^l + — |1- 



and 



e(fc,w)~l + — 



1 - 



— k 2 i 



u) 2 ^ — k 2 vp — ij 



(A7) 



(A8) 



These are precisely the same DFs obtained previously in Refs. [23, [28|, l23| which used the same small k,u>- 
approximation limits of the more general expressions (|A1[) and (|A3[) . Note that the DFs (|A7[) and (|A8p are the 
quasiclassical limits of the more general expressions (|A1[) and (|A3|) , respectively. Therefore they can be alternatively 
derived from a classical kinetic equation within RTA with the Fermi distribution function as an unperturbed state. 



APPENDIX B: EVALUATION OF THE PARAMETERS hi AND h 2 



In this Appendix we give detail derivation of the parameters hi and hi which contributes to the high-velocity SP 
of an interacting 2D electron gas, Eqs. and ([IS]) , respectively. First we write Eq. (jTTjl in another but equivalent 



form, hi = p(s)\ s 



where 



pO) 



7T 

8a 

^2 



A (u) du — 

OO PS 

zdz I L (z, u) udu - 



(Bl) 



zdz I L [z, u) udu 



Here L(z,u) — Im[— l/e(z, u)] is the energy loss function. For derivation of Eq. (jBlj) the Bethe sum rule (the first 
frequency moment of the energy loss function) in variables z and u has been used (see, e.g., Refs. 0,0,11]) 



Im 



-1 , *x z 

— -udu = 

e(z,u) 4z 



7T 

Aaz 



Assuming that the upper cutoff s is large enough, s> 1, Eq. (|B1[) can be written in explicit form 



p{s) 



8a 



s+l 



zdz 



L (z, u) udu - 



z-l 



zdz 



s-l 



z+1 



L (z, u) udu 



(B2) 



(B3) 



16a 2 



g(z) z 3 g 2 (z)u r (z)dz 

\9{»)\ \<Pr(z)\ 
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Here u r (z) is the solution of the dispersion equation for an interacting DEG, Eq. (|16p . and we have introduced a lower 
cutoff parameter z r (s) which z, r (s) — > at s — > oo. Also we have introduced the function <p r (z) which is given by 



d 

§r (z) = Tj-fl (Z,U) 



u r (z) — z 



u r (z) + z 



,(») 



\J[u r (z) - zf - 1 y [u r (z) + zf - 1 



(B4) 



The last term in Eq. (|B3|) is the contribution of the collective excitations and hence the function fi(z, u) in Eq. (|B4|) 
is defined in the domain < z < u — l.(or u > z + 1). Without lose of the generality we chose as a lower cutoff (i.e. 
z r (s)) a function which is inverse to u r (z). Using Eq. p^|) it is straightforward to calculate the asymptotic behavior 
of this function at large u. It behaves as 



z r (u) 



2au 2 



it 2 u 4 



(B5) 



where the expansion coefficients are given by 

At 



3 1 

4 -2 X ( r «)- 



^2 = -[l-x(r s 



x(r s ) 
16 



[3-x(r s )], 



(B6) 
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A, = ~T 



1 



64 4a 2 



x(r s ) 
16 



29 

x 2 (r s ) - 2x 2 (r s ) + 9x(r s ) - — 



with x(r s ) = G (r s ) /[aG 12 (r s )] and x (r s ) = G 22 (r s ) /[aG 12 (r s )}. 

Since at small z the functions u r (z) and 4> r (z) behave as u r (z) ~ (2az)- 1 / 2 and 



<}> r (z) ~ 4aV2^z 5 / 2 1 1 + ^ [1 + 2x (r s )] + O (z 2 ) j 



(B7) 



respectively, at s — > oo the plasmon contribution in Eq. (|B3|) vanishes as ~ z r (s) ~ s -2 — > 0. For calculation of the 
first two terms in Eq. (|B3[) (single-particle contributions) we first make a substitution of the integration variables, 
z — » z + s and u — ► it + s. At s — > oo the remaining expression behaves as p (s) ~ —1/ (27rs) — > 0. Thus at s — > oo, 
p (s) — > and 



hi = p(s)\ ( 



= 0. 



(B8) 



For calculation of the coefficient h 2 it is imperative to evaluate the third moment of the energy loss function. In 
2D and in general case this has been done in Ref. @. In the present context of an interacting electron gas with DF 
([T3")) and (fi"4"l) this moment is given by 



Im- 



-1 



du =£-<j-l-[l- G(z)} + 1 + 2; 



e{z 1 u) 8a [az 2 L WJ ' 2z 

As we have done above we introduce now a new function through the relation h 2 = 3 (s) where 

s s 3 1 
4G^~ 3G^~ 2 

and C3 is given by Eq. (J39J) . Employing the relation (|B9j) for the third frequency moment we obtain 

1 1 

2 



1 f s 

In 3 (s) = — / A (u) u 2 du — In ; 
C3 Jo 



In 3 (s) 



3 [1 - xi (r»)] 
1 

_ 2[l-x!(r s )] 



[In (2a) +3xi (r s )lnx 2 (r s )] 



(B9) 



(BIO) 



(Bll) 



E/(s), 



where Xi (r s ) = Go (r s ) /[3G22 (?"s)]> X2 (r s ) = G12 (r s ) /[2G 22 (r s )] are new density-dependent parameters. The 
function U(s) is evaluated in the similar way as we have done above. In particular, neglecting the contribution of 
plasmons which is again vanishingly small at s — > 00 this function becomes 



U(s) 



16a 2 
~3tT 



zdz 



z+l 



L (z, u) u 3 du - 



f S+1 f S ITS 

/ zdz / L (z,u) u 3 du — — 
Js Jz-l 8a 



(B12) 
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Now only the single-particle excitations contributes to Eq. (|B12[) . Again by making the changes of the integration 
variables, z — ► z + s and u — ► u + s, at s — ► oo we have found that the function U(s) behaves as 



U{s) = --k 1 (r s ) + - + 0(s- 2 ). 
3 os 

Finally, substituting Eq. (|B13|) into Eq. (|B11[) and taking the limit s — > oo we arrive at 



(B13) 



lnh 2 (r s ) = 



1 



3[l-xi(r a )] 



In 



2^2 



+ (r s ) [1 + ln^ 2 (r s )] - 2 



(B14) 



The transition to the limit of non-interacting 2D electron gas is performed by taking the limit x\ — > in Eq. (|B14[) 
which yield 



ln/i 2 (r s ) 



In 



'2%/2 N 



(B15) 
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